Off-shell effects and consistency of many-body treatments of dense matter 
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Effective field theory requires all observables to be independent of the representation used for the 
quantum field operators. It means that off-shell properties of the interactions should not lead to any 
observable effects. We analyse this issue in the context of many-body approaches to nuclear matter, 
where it should be possible to shift the contributions of lowest order in purely off-shell two-body 
interactions into three-body forces. We show that none of the commonly used truncations of the two- 
body scattering amplitude such as the ladder, Brueckner-Hartree-Fock or parquet approximations 
respect this requirement. 
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Nuclear forces cannot yet be derived from QCD, and so 
we must rely on phenomenological nucleon-nucleon po- 
tentials There is a variety of these potentials, all 
of which give similarly good fits to the available low- 
energy two-body scattering data but which have differ- 
ent off-shell behaviours. In recent years it has become 
clear that a consistent description of nuclear systems re- 
quires three-body forces for both systems consisting of 
a small number of nucleons ||, as well as for nuclear 
matter j|. Part of the role of the three-body force is to 
compensate for the different off-shell behaviours of the 
two-body forces. This is as expected because the physics 
should not depend on the off-shell behaviour of the inter- 
actions. Exact calculations with phase-equivalent two- 
body forces and their corresponding three- (and higher-) 
body forces should therefore give identical results. How- 
ever this leaves open the question of whether the approx- 
imate many-body techniques used in actual calculations 
respect this property. 

The most elegant way to analyse this problem is to 
use effective field theory (EFT), which has become an in- 
creasingly popular tool in modern nuclear physics. This 
approach makes use of the fact that the dynamics at low 
energy is only weakly dependent on the high-energy de- 
grees of freedom, and so a detailed knowledge of the inter- 
action at short distances is not required. The low-energy 
physics can then be described using a local effective La- 
grangian jjj. The physical amplitudes can be obtained 
from this in the form of expansions in powers of the low- 
energy scales involved. 

Application of the EFT approach to the nucleon- 
nucleon (NN) system turns out to be rather complicated 
due to the large s-wave scattering length, which gives rise 
to an additional small scale in the problem. If the scat- 
tering length were similar to the range of the NN inter- 
action (usually referred to as "natural" ) then it would be 
possible to use a perturbative treatment. Instead, the un- 
naturally large scattering length means that the leading 
terms in the chiral expansion of the NN potential should 



be iterated to all orders, as shown by Weinberg (D an d 
van Kolck ||. The remaining terms in the potential can 
be treated as perturbations, organised according to the 
power counting elucidated by Kaplan, Savage and Wise 
(KSW) 0. This can also be thought of as an expan- 
sion around a renormalization-group fixed point which 
corresponds to a bound state at threshold || . 

Since the off-shell properties are determined by the 
representation chosen for the field operator in EFT's, 
independence of physics on off-shell properties is called 
"reparametrization invariance" . In formal field theory 
the corresponding result is known as the equivalence the- 
orem H . This issue has been considered in several recent 
papers juj, @. (Further refer ences can be found in 
Ref. pflT) Furnstahl et al. |llj have demonstrated that 
for a model consisting of a dilute Fermi system with a 
natural two-body scattering length three-body counter- 
terms can indeed cancel the effects of off-shell part of the 
two-body interaction. 

Here we explore the extension of these ideas to dense 
systems with an unnatural scattering length, such as nu- 
clear matter. We find that the effects of off-shell two- 
body interactions can indeed be cancelled by three-body 
forces. This is unsurprising given the general nature of 
the equivalence theorem. Of more importance are the 
sets of diagrams that need to be included to obtain this 
cancellation, and the implications for commonly used 
many-body truncation schemes. We shall use our analy- 
sis to show inconsistency of various many-body methods. 

At low enough energies we can work with EFT's for the 
NN system where the interaction can be treated as purely 
short-ranged. In principle these should be improved by 
by including one-pion exchange explicitly, but there is 
still some debate about the best way to do this. (See, for 
example, Ref. [[l3] and references therein.) We consider 
first an EFT which leads to a purely energy-dependent 
NN potential. To next-to-leading order (NLO) in the 
small scale this spin- and isospin-independent potential 



has the form 



V 1 = C + C' 2 p 2 



(1) 



where p 2 /M is the relative kinetic energy of the nucleons. 
The LO coupling constant Cq is of order (where Q 
is a generic low-energy scale) in KSW counting and so 
it should be treated nonperturbatively 0. The NLO 
coupling C 2 is proportional to the effective range, and 
is thus of order Q° and can be treated perturbatively in 
this counting scheme. Nonetheless, to show an example 
of the principles discussed above, we solve the Lippmann- 
Schwinger (LS) equation with this potential to all orders 



in C 2 , and find the vacuum T-matrix 



C +P 2 C 2 



1 + H(Co+p 2 ^)(zp + m) 



(2) 



Here we have used a subtractive renormalization pro- 
cedure [frl (All coupling constants here should be 
understood as renormalized ones which depend on p to 
ensure that the scattering amplitude is /^-independent.) 

More generally, the effective Lagrangian can also in- 
clude interactions with space derivatives of the nucleon 
fields and this leads to a potential that depends on mo- 
mentum as well as energy. The most general NLO po- 
tential has the form 



V 2 = Co + C' 2 p 2 



\c 2 (k 2 



k' 2 -2p 2 ), 



(3) 



where k and k' denote the initial and final relative mo- 
menta of the nucleons. The coupling C 2 describes a 
purely off-shell interaction. Solving the LS equation we 
get 



1 + 



1 
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k 2 )( P 2 



C 2 (k 2 +k' 2 -2p 2 



k' 2 ){ip + v) 



,(4) 



where T\ is given by Eq. (g). From this we can see that 
the two T-matrices coincide on-shell (k 2 = k' 2 = p 2 ) and 
so the scattering observables are indeed independent of 
the off-shell behaviour of the potential as required by the 
equivalence theorem. 

The situation becomes much less trivial in the pres- 
ence of the nuclear medium. An in-medium T-matrix 
Jill can be obtained by solving the Feynman-Galitskii 
(FG) equation, 



T m = V + VG F T r ' 



(5) 



where G F denotes the two-nucleon propagator which con- 
tains both particle-particle {pp) and hole-hole (hh) states. 
This T-matrix can be thought of as an extension of the 
more familiar G-matrix |L6|, [llj to include hh as well as 
pp ladders. It is convenient to represent the FG equa- 
tion graphically in terms of the Hugenholtz diagrams of 





FIG. 1: Hugenholtz diagrams representing the equation for 
(a) the T matrix and (b) the dressed propagator. The open 
circle denotes the LO potential Co, and the solid circle the T 
matrix. 



Fig. ga). These are versions of Feynman diagrams which 
explicitly incorporate antisymmetry of the interactions. 
Internal lines represent Feynman propagators which de- 
scribe both particles and holes. The arrows represent the 
flow of quantum numbers such as baryon number. Each 
topologically distinct diagram should be multiplied by a 
symmetry factor to take account of the number of ways 
it can be constructed from the antisymmetric vertices. 
More details of these diagrams and the rules for evaluat- 

[nil. 



ing them can be found in the textbooks 16 



The solution of the FG equation is rather straightfor- 
ward in the case of zero total momentum of the nucleons 
For the potential V\ it takes the form 

1 



M 
4ir' 2 



(6) 



where pf is the Fermi momentum. In the same way we 
can solve the FG equation for the potential V 2 . We shall 
assume that the C 2 term can be treated as a perturba- 
tion. For simplicity we omit the energy-dependent term 
C 2 from now on. Although this term makes a physically 
important contribution to the energy of the two-particle 
amplitude, it does not take part in the cancellation of 
off-shell effects which is of interest here. To first order in 
C'2 the in-medium T-matrix can be written 



rptn rprr, 

1 2 — 1 1 



rrg(2 P 2 



k' 2 ) 



(7) 

If we now evaluate T™ at the on-shell point, we see that 
it does not agree with Tf since the last term does not 
vanish. This indicates that calculations of nuclear matter 
based on the in-medium T-matrix (or similarly the G- 
matrix) do not satisfy the requirement of reparametriza- 
tion invariance. Alternatively, in more traditional nu- 
clear physics language, results for in-medium observables 
depend on the off-shell behaviour assumed for the NN po- 
tential. Such a dependence is unphysical and should not 
be present. A clue to how the dependence may be re- 
moved comes from the form of the final term in Eq. (fif), 
which is proportional to the density. Its structure is thus 
similar to that arising from a three-body contact inter- 
action. This suggests that it may be possible to trade off 
the off-shell dependence against a three-body force. As 
shown below, this can be done, provided our approach 
includes more than just ladder diagrams. 
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FIG. 2: Hugenholtz diagrams for the ground state energy at 
first order order in C2 (the open triangle). 




FIG. 3: Hugenholtz diagrams for the ground state energy at 
first order in the three-body force Do (the open square). 



Before exploring what additional physics is needed to 
remove the off-shell dependence, we should note that 
there is no clear separation of scales in strongly inter- 
acting, dense matter. This is an unsolved problem for 
the application of EFT's: no power counting has been 
found which leads to a consistent expansion. One should 
really solve the many-body theory for Co exactly, by con- 
structing the full in- medium NN vertex, T. Nonetheless 
simpler approximations are commonly used in nuclear 
physics, typically replacing the full NN vertex by a G- or 
T-matrix. Including ph rings as well as pp and hh ladders 
leads to the parquet approximation |l6|, |l8). We exam- 
ine here the consistency of these approximations with 
reparametrization invariance. 

At LO in C 2 the contributions to the ground-state en- 
ergy of matter are shown in Fig. |^, where the solid dot de- 
notes an in-medium NN vertex. If Co were weak enough 
we could expand these diagrams perturbatively to get a 
contribution of order MC0C2. The resulting diagrams 
have an identical structure to those in Fig. |] below, ex- 
cept that none of the propagators are dressed. As shown 
in Ref. jyj, they can be exactly cancelled against the 
LO contribution of a contact three-body interaction with 
strength Dq = Y1MCqG 2 . This is as required by the 
equivalence theorem, since the off-shell term and three- 
body force with this strength are both generated from 
a Lagrangian which contains neither by the same field 
redefinition. The details are given in Ref. pj . For defi- 
niteness we repeat the relevant Feynman rules here: the 
two-body vertices, represented by an open circle and a 
triangle respectively, are 

- iC S 2 and iC 2 (A l + A e + Aj + A j ,)S 2 , (8) 

where A, = Mpf — (pi) 2 /2, pi being the four- momentum 
of the ith nucleon, and the spin-isospin structure is given 



by S 2 = Su'Sj 



The three-body vertex (an open 



square) is 



FIG. 4: The contributions from Fig. Q proportional to the in- 
tegral Iq. Diagrams which can be obtained from those shown 
by simply reversing all the arrows are not shown separately. 





FIG. 5: (a-c) Contributions proportional to the integral Ii 
obtained from Fig. || in the BHF approximation, (d) Ex- 
tra contribution from the parquet approximation. Diagrams 
which can be obtained from these by simply reversing all the 
arrows are not shown separately. 



When the leading two-body vertex is resummed we get 
an effective vertex T(pi , pj ; p^ , pj/)S 2 which is denoted by 



a filled circle 20 



Treating Co nonperturbatively, the same three-body 
force gives rise to the diagrams shown in Fig. ||. Each 
of these diagrams gives a contribution equal to a distinct 
integral multiplied by Dq and a degeneracy factor. The 
detailed forms of these integrals, which we denote by Iq, 
ii, I 2 , ■ ■ ■ are not needed here. However we can use these 
integrals to classify the structures which arise from the 
diagrams of Fig. |[ To evaluate them, we note that the 
off-shell vertex can be written as a sum of four pieces, 
each of which can cancel a bare propagator on one "leg" : 



G (q)(Mq - q 2 /2)C 2 = iMC 2 , 



(10) 



iD 



Su'(5jj'hk> - 6jk'S k f) + cyclic(i', j',k') . (9) 



where Gg(q) is the bare single-particle propagator. The 
diagrams in Fig. |^ give rise to many different contribu- 
tions, which can be identified by iterating the equations 
for the in-medium NN vertex and dressed propagator to 
pull out a bare propagator ending on a lowest-order ver- 
tex Co on any of the lines in the original diagrams. When 
the bare propagator is cancelled against the off-shell ver- 
tex, as in Eq. (|l0|), the result is one of the integrals /„ 
multiplied by MCqC 2 and a numerical factor. Thus we 
can examine the cancellation of off-shell dependence for 
each integral in turn. 

We consider first the Brueckner-Hartree-Fock approxi- 
mation (BHF) |ll| [ijj , in which propagators are dressed 
and the in-medium NN vertex is obtained by iterating the 
potential in the pp and hh channels, as shown in Fig. [I]. 
The contributions proportional to Iq from Figs. |^(a) and 
(b) are shown in Fig. |^. Except for the dressing of the 
propagators, these have the same structures as the per- 
turbative diagrams considered in Ref. |Tl"] ] and they can 
be shown to cancel with Fig. |3|(a) in the same way. 
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FIG. 6: (a) Contribution from Fig. y(b) proportional to the 
integral I2 in the BHF approximation, (b-c) Extra contribu- 
tion in the parquet approximation, (d) Diagram containing a 
non-parquet contribution. Diagrams which can be obtained 
from those shown by simply reversing all the arrows are not 
shown separately. 



Fig. ^(a) gives one other contribution, shown in 
Fig. g|(a), which is proportional to I\. In the ladder ap- 
proximation to the NN vertex, Fig. ||](b) also gives con- 
tributions proportional to ii, shown in Fig. |^(b-c). The 
sum of Figs. |(a-c) is -2g(g-l)(2g-3)MC C 2 Ii, where 
g is the spin-isospin degeneracy factor (g = 4 for sym- 
metric nuclear matter). In contrast, the three-body force 
gives g(g - l)(g - 2)D h/2 from Fig. |(b). We see that 
the degeneracy coefficients do not agree and cancellation 
does not occur. For example, the off-shell dependence 
is nonzero for neutron matter (g = 2) where the Pauli 
principle forbids a contact three-body force. 

There is one other structure proportional to ii, 
Fig. ||(d). However this cannot be generated from the di- 
agrams of Fig. if the potential is iterated in the pp and 
hh channels only; it requires iteration in the ph channel as 
well. When this contribution is included, the degeneracy 
factors agree and the off-shell dependence proportional 
to 1\ is indeed cancelled by the three-body force with 
Do = 12MC C 2 . 

The crucial point to note is that the cancellation re- 
quires diagrams which can only be obtained by iterating 
the two-body potential in the ph channel. These are not 
contained in the ladder or BHF approximations and so 
any approach based on a G- or T-matrix cannot satisfy 
the equivalence theorem. Observables calculated in these 
approaches will have an unphysical off-shell dependence 
which cannot be absorbed into a three-body force. 

The need for diagrams with iteration in the ph chan- 
nels suggests that one should try a more complete ap- 
proach. One such, which treats all two-body channels in 
a symmetric way, is the parquet approximation 
If we interpret the solid circles in Fig. || as parquet NN 
vertices constructed from Co, then all of the contribu- 
tions in Fig. H can be generated by iterating the parquet 
equations. (Note that the parquet self energy can still 
be expressed in the form of Fig. |l|(b) p9[ , and so the 
discussion of 1$ above is unchanged.) 

Turning now to terms proportional to I2, which ought 
to cancel with the three-body graph Fig. |^(c), we find 
one self-energy contribution, Fig. ||(a). This would be 
present even in the BHF approximation for the T-matrix 
but, not unexpectedly, this does not provide the cancel- 



lation. In the parquet approximation, Fig. ||(b) gives 
additional contributions, shown in Fig. ^|(b,c), and part 
of (d). Only when the non-parquet contributions are in- 
cluded (see Table 3 of Ref. [jl8| ) , i.e., with a full set of 
diagrams, do we match the result from Fig. ||(c). Thus 
we conclude that the parquet approximation also violates 
rcparamctrization invariance! 

It may well be possible to include the necessary struc- 
tures by extending the parquet approximation along the 
lines discussed in Ref. | 16fl , starting from a basic vertex 
which is a sum of diagrams which are two-particle irre- 
ducible in all channels. However if, as suggested there, 
these structures are simply added in perturbatively, they 
will not generate the full in-medium vertices needed for 
the diagrams of Fig. |[ 

In summary, our results demonstrate that the require- 
ment of reparametrization invariance, which would re- 
quire the effects of off-shell dependence of the two-body 
interaction to be cancelled by a three-body force, are 
not satisfied by any of the commonly used truncations of 
the two-body scattering amplitude such as ladder, BHF 
or parquet approximations. The violations show up in 
structures with higher numbers of insertions of the in- 
medium NN vertices for the more sophisticated trun- 
cations, but as the interaction is strong this does not 
provide a consistent expansion scheme. Finding such 
a scheme remains essential for practical applications of 
EFT's to dense, strongly interacting matter. 
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